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Abstract 



We have recently derived two-fluid hydrodynamic equations for a trapped 
weakly-interacting Bose gas. In this paper, we use these equations to discuss 
first and second sound in a uniform Bose gas. These results are shown to 
agree with the predictions of the usual two-fluid equations of Landau when 
the thermodynamic functions are evaluated for a weakly-interacting gas. In a 
uniform gas, second sound mainly corresponds to an oscillation of the super- 
fluid (the condensate) and is the low frequency continuation of the Goldstone- 
Bogoliubov symmetry-breaking mode. 

PACS numbers: 03.75.F, 67.40.Db 
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I. INTRODUCTION 



One of the most spectacular features |TJ exhibited by superfluid 4 He is the existence of 
two hydrodynamic sound modes, first and second sound. As first pointed out by Tisza 0, 
the motion of a Bose condensate as a separate degree of freedom results in a two fluid hy- 
drodynamics describing the superfluid and normal fluid components ||. In recent work, the 
authors Q gave a microscopic derivation of the two-fluid hydrodynamic equations of motion 
for a trapped weakly-interacting Bose-condensed gas. In contrast to a Bose-condensed liq- 
uid like superfluid 4 He, the superfluid in a gas corresponds directly to the condensate atoms 
and the normal fluid corresponds to the non-condensate (or excited) atoms. In the present 
paper, we use these two-fluid equations to discuss the first and second sound modes of a uni- 
form Bose-condensed gas. We find that at temperatures close to Tbec, first (second) sound 
mainly corresponds to an oscillation of the non-condensate (condensate) atoms. We also 
confirm [^J|| that it is the second sound mode in a uniform gas which is the low frequency 
hydrodynamic analogue of the collisionless Bogoliubov- Popov Goldstone mode 0. 

These results for a uniform gas are of interest for comparison with the hydrodynamic 
oscillations of the condensate and non-condensate in a non-uniform trapped Bose gas ||||. 
They also may be of direct interest in connection with recent studies at MIT of the 
propagation of pulses along the z-axis of a cigar-shaped trap. The axial trap spring constant 
is so small that the condensate along the z-axis can be treated as effectively uniform (to a 
first approximation) in such propagation studies. 

We recall that Ref. |4| (ZGN) is based on: (a) a time- dependent Hartree-Fock-Popov 
equation of motion for the condensate wavefunction $(r,t); and (b) a set of hydrodynamic 
equations for the fluctuations of the thermal cloud (non-condensate) based on a kinetic equa- 
tion which includes the effect on the atoms of the time-dependent self-consistent Hartree- 
Fock field. The analysis of ZGN uses the local equilibrium solution of the kinetic equation 
and thus does not include any hydrodynamic damping, such as Kirkpatrick and Dorfman [H^ 



consider. However it should be emphasized that a local equilibrium description is crucially 
dependent on collisions between the atoms and thus the hydrodynamic equations are only 
valid for low frequency phenomena (w < 1/ r c , where r c is the mean time between collisions 
of atoms in the thermal cloud). 

In Section |TT], we solve the linearized hydrodynamic two-fluid equations for the coupled 
superfluid and normal fluid velocity fluctuations derived in ZGN. We exhibit the first and 
second sound normal modes valid at intermediate temperatures, defined as the temperature 
regime below Tbec where the interaction energy of an atom is much less than the thermal 
kinetic energy (i.e., gn <C k^T; here, n is the gas density and g = Anatf/m is the 
interaction parameter). The analysis of ZGN is built on a mean-field approximation for the 
equilibrium properties. As discussed in Section III , this simple theory is not valid close to 
the superfluid transition, where it gives rise to spurious discontinuities in the condensate 
density. In Section [TV], we discuss the relation between our two-fluid equations written in 
terms of velocity fluctuations and the standard Landau formulation given in terms of density 
and entropy fluctuations [|1],[3] . All previous discussions |TT||7l,p|1 of hydrodynamic modes in 
a dilute Bose gas have used the latter formulation. 



2 



II. COUPLED EQUATIONS FOR SUPERFLUID AND NORMAL FLUID 

VELOCITIES 



When there is no trapping potential, the non-condensate density h and condensate 
density n c0 do not depend on position. In this case, one can reduce the linearized 2-fluid 
equations given by Eqs. (12), (15) and (16) of ZGN to two coupled equations for the normal 
and superffuid local velocities 
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We emphasize that these equations are only valid at finite temperatures such that grio <C 
ksT. In deriving these equations, we have assumed that the contribution from the first term 
of Eq.(13) of ZGN is negligible in the long- wavelength limit of interest. These equations 
can be solved to give the low frequency hydrodynamic normal modes of a uniform Bose- 
condensed gas, as will be discussed. We defer discussion of the equilibrium quantities (n c0 , 
ho and the kinetic contribution to the pressure Po) which appear in (|la|) and ([TtJ) to Section 
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The zeros of the secular determinant of this coupled set of equations give two phonon 



solutions uj± 



u±k 2 , where the velocities are the solution of 
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Expanding to second order in the explicit dependence on g, the sound velocities are given 
by 
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where e = Agrio/^^ - <C 1 is the expansion parameter. We note (see Section III) that the 

O Tlo 

ratio = fcgT ^ 2 /^ ) depends weakly on g. The uo + mode in (pJ) clearly corresponds to 
n o g3/2{Zo) 1 — 

first sound. Using uj\ = u 2 ,k 2 in (0), one finds to leading order in g that 
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^ - - » 1 ■ (5) 

That is to say, the uj + first sound mode corresponds to an in-phase oscillation in which 
the non-condensate velocity amplitude is much larger than that of the condensate. The uj- 
mode in (|4"15 ) is the second sound mode. Using lo 2 _ = u 2 _k 2 in (^|), one finds to leading order 
in g that 

<Ps 2 ho 

- — ~ > 1 . (6) 

(Pn e n c0 

Thus at finite temperatures where (|la|) and (|Tb"l) are valid, second sound in a uniform weakly- 
interacting gas is seen to be an out- of -phase oscillation, in which the condensate velocity 
amplitude is much larger than that of the non-condensate (a similar result was obtained 
many years ago in Ref. 



III. EQUILIBRIUM PROPERTIES IN THE POPOV APPROXIMATION 

The ZGN derivation of the coupled hydrodynamic equations for the two velocity fields 
given in ( |ia| ) and ( [TBI ) is built on a self-consistent Hartree-Fock description of the equilibrium 



properties. One of the earliest discussions of this mean-field theory was given by Popov 
and it has become the standard approximation in recent studies of trapped Bose gases. 
Referring to ZGN, we recall that the equilibrium equation for the condensate yields the 
equilibrium chemical potential 

Ho = 2gn + gn c0 . (7) 

This parameter enters in the determination of the equilibrium excited-atom density given 
by 

n (T,n Q ) = -j^g 3/2 {z Q ) , (8) 



Zq — e /30o-2gn ) _ e -Pgn c0 ^ 



where (n = n c0 + n ) 



is the equilibrium fugacity and A = y2irh 2 /mksT is the thermal de Broglie wavelength. 
The associated excited-atom kinetic pressure is 

> n o) — ^ 3 5'5/2(^o) • (10) 



We note that these results are equivalent to the simple "toy model" studied in Ref. [JT2 . 

Eqs.(H) and (|[) must be solved self-consistently to determine n c0 and h for a given 
total density n . Condensation occurs when the density reaches the critical density n cr = 
^3/ 2 (l)/A 3 . For n < n cr , the condensate density is zero and (|J) with h = n determines the 
equilibrium fugacity. In Fig. 1, we show the equilibrium densities as a function of volume for 
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a fixed temperature. The parameter 7 cr = (3gn cr is used to characterize the strength of the 
interaction. We see that the present level of approximation leads to a discontinuous change 
in the densities at the transition point [D|. Moreover, below the critical volume v cr = l/n cr , 



ho decreases as a result of the interactions with the condensate, in contrast to the ideal gas 
behaviour which has the non-condensate maintaining a constant density of n cr . Fig. 2 gives 
the total pressure defined as [§J 

P = P + ^g(n 2 + 2n h - hi) , (11) 

normalized by the critical pressure P cr = g 5 / 2 {l)/ (3A 3 of the ideal gas. The second term 
in (|TT|) is the explicit interaction contribution, but it should be noted that Pq also depends 
on interactions as a result of its dependence on zq. The discontinuous behaviour of the 
non-condensate density leads to an analogous discontinuity in the pressure In Figs. 



3 and 4 we show the corresponding behaviour as a function of T. It is of interest to note 
that for a trapped Bose gas, the use of these equilibrium properties in the Thomas-Fermi 
approximation leads to a similar discontinuous behaviour of the equilibrium condensate 
density, but now as a function of the radial distance from the center of the trap [14]. 



It is clear that the properties of the weakly-interacting gas are nonanalytic functions of 
the interaction strength g at the transition point within the mean-field Popov approximation 
described by (|7|) (JT0|) . However, one should not take these features in the BEC critical region 



seriously. The simple mean-field Popov approximation for interactions is well known ||15|116 
not to be valid very close to the transition and the predicted discontinuities exhibited in 
Figs. 1-4 (characteristic of a first-order transition) are indicative of the limitations of the 
present simple theory. A correct treatment of this region would require a renormalization 
group (RG) analysis JT7| which is outside the scope of the present paper. 

For later purposes, we note that the kinetic pressure P in (|T^) can be calculated by 
expanding the fugacity as z ~ 1— f3gn c0 + - ■ -, which yields (using the identity zdg n (z)/dz = 
g n -i(z)) 

h - Per ~ gn c0 n CT , (12) 

where P cr and n cr are the critical pressure and density of the ideal Bose gas introduced earlier. 
However, a similar perturbative expansion of the non-condensate density h in (|§|) is not 
possible since the derivative of gz/2{z) diverges at z = 1. Indeed, it is this non-perturbative 
dependence on g which leads to the discontinuities shown in Figs. 1-4. 



IV. RELATION TO STANDARD TWO-FLUID EQUATIONS 

First and second sound in a uniform Bose-condensed gas have been previously discussed 
in the literature [^,f7j|lT| . These earlier treatments start with the usual two- fluid equations 
of Landau 0. We recall that these linearized equations are (see ch.7 of Ref. |L|) 
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where 



8n(r, t) = 5h(r, t) + 5n c (r, t) 
5j(r, t) = n 5v N + n c0 5v s ■ 



(14) 



P and s are the pressure and entropy density, respectively. ZGN proved that the two-fluid 
equations which lead to ( |TaD and ([l|]) are in fact equivalent to the two-fluid equations in 
(O) when the thermodynamic functions in the latter are evaluated for the present model of 
a weakly-interacting Bose gas. Using the thermodynamic relation ||], n$5p = <5P — Sq5T, to 
eliminate the chemical potential, and defining the entropy per unit mass by s = s/mn = s/p, 
one can reduce the equations in ( |13|) to [l],||] 
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Solving this closed set of equations in terms of the variables Sn and Ss, one finds two normal 
mode solutions u 2 = u 2 k 2 , where u 2 is given by the solution of the quadratic equation [EJ 
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In this equation, c v is the specific heat per unit mass and derivatives of the pressure have 
been expressed in terms of the independent thermodynamic variables T and p. Although 
not immediately apparent, the coefficients in (|i~6"D are in fact consistent with those appearing 
in©- 

The problem is thus reduced to evaluation of the various equilibrium thermodynamic 
functions and derivatives which appear in (|TE|). For the entropy per unit mass we have the 
expression [|| 
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from which we obtain 
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From the equation of state (IT), we find that 
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and 
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These quantities have been calculated previously in the limit that the interaction parameter 
g is regarded as small 0,0,0. In this situation, P in ( [TP] ) is approximated by flT2"|). An 
additional approximation is typically made whereby ho is simply replaced by the ideal gas 
expression n cr , in which case n c o = no — n cr . To the same level of approximation, one finds 
|s = and ^ = 3n cr /2T. With these replacements, we also note that the expressions 

for the pressure and the entropy and energy densities given in ZGN reduce precisely to those 
of Refs. @ and [0 . 

Using these results to calculate the thermodynamic quantities in flT6|) , the first and second 
sound velocities are found (after some algebra) to be given by 
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keeping terms to first order in g. The leading order terms in ( [2 lap and ( pib| ) were obtained 
from (|i~6|) by this method by Popov (see the last paragraph of Ref. [0]) as well as by Lee 
and Yang [XT]. Precisely the same results follow from (^) to first order in g when (|T2| ) is 
again used for the kinetic pressure Po and h is replaced by n cr . However, the results given 
by (^) are more general than those in (|2"T|), which only keep the leading order corrections to 
the properties of a non-interacting gas. As we discussed above, the analysis leading to (|2T|) 
ignores any interaction-correction to the non-condensate density h, which as can be seen 
from Fig. 1, becomes significant as the density increases beyond n cr . 

As we emphasized in the beginning of Section ||, the analysis of ZGN assumes that 
gn Q <C k B T and thus the results are not really valid at low temperatures. To discuss the 
low temperature region would require a generalization of our work which is based on a 
quasiparticle spectrum exhibiting phonon-like behavior at long wavelengths (such a kinetic 
equation has been derived in Ref. [K|). The pioneering work of Lee and Yang did 
include an analysis of both the low temperature and high temperature regions. At low 
temperatures, they found that the first and second sound modes avoid becoming degenerate 
by hybridizing and an interchange of the physical meaning of these two modes occurs as a 
result of this hybridization. While the sound velocities given by (^|) are not really valid at 
low temperatures, Fig. 5 shows that our results do lead to this expected hybridization of 
first and second sound in a dilute gas. 



V. CONCLUDING REMARKS 



Recently two-fluid hydrodynamic equations were derived [f| for a trapped, weakly- 
interacting Bose gas. These are given in terms of coupled equations for the superfluid 
and normal fluid velocity fluctuations. In order to obtain more physical insight into these 
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hydrodynamic equations, we have given in the present paper a detailed analysis for a uni- 
form Bose gas. In this case, it has been proven that the hydrodynamic equations of Ref. 
H are formally equivalent to the usual Landau two-fluid equations. As the present paper 
shows, this formal equivalence is somewhat hidden in explicit calculations of the first and 
second sound velocities. However, as discussed in Section [TV], our results do reduce (to first 
order in the interaction g) to those found in earlier studies p],|7||5| based on the Landau 
formulation. 

In superfluid 4 He, one evaluates the equilibrium thermodynamic parameters in ( |16D us- 
ing the phonon-roton excitation spectrum. As is well known [0,0, in superfluid 4 He, first 
sound corresponds to an in-phase oscillation in which vjy = V5. In contrast, second sound 
corresponds to an out-of-phase oscillation in which p ri yN = —ps^s- The difference between 
second sound in a dilute gas at finite temperatures and in a liquid is a result of the dom- 
inance of the kinetic energy over the interaction energy for atoms in a gas. In both cases, 
however, we note that the second sound frequency goes to zero (becomes soft) at the super- 
fluid transition. The mode does not exist above Tbec- Moreover (|4b|) shows that second 
sound crucially depends on the interaction g. It would be absent if we had set g = in ( [Taj) 
and dlb]). 

As we have noted, second sound in a dilute gas largely involves an oscillation of the 
condensate atoms (superfluid density) and is a soft mode which vanishes in the normal 
phase. We recall that at finite temperatures 0, the generalization of the T = Bogoliubov 
phonon gives a velocity formally identical to the first term in (|4b|) . Thus we conclude that 
in a weakly interacting Bose-condensed gas at finite temperatures, second sound is the low 
frequency (hydrodynamic regime) continuation of the high frequency (collisionless or mean- 
field regime) Bogoliubov-Goldstone mode. This was first suggested in Refs. 0-0. The 
situation is quite different in superfluid 4 He, where the collisionless phonon spectrum is the 
continuation of hydrodynamic first sound and there is no high-frequency analogue of the 
second sound branch. 
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FIGURE CAPTIONS 



Fig.l: Density vs. volume per particle for a fixed temperature T. The chain curve corresponds 
to the non-condensate, the solid curve to the condensate. 7 cr is the value of gn /kBT 
at the critical density n cr = #3/2(1) /A 3 . 

Fig.2: Pressure isotherms: the solid line is the total pressure according to flTTp, the chain 



curve is P and the dashed curve corresponds to the usual approximation [|T8],|T9| P 

P cr + \g{n 2 + n 2 cr ). 

Fig.3: As in Fig. 1, but as a function of T for a fixed density n . Here, 7 cr = gno/ksTBEC- 

Fig.4: Normalized pressure as a function of T for a fixed density n . The solid curve corre- 
sponds to ( PD and the chain curve is P . The dashed curve below T = T BEC is the 
ideal gas result P /P cr = (T/T BEC ) 5/2 . 

Fig. 5: Squares of the first and second sound velocities (normalized by the first sound velocity 
of the ideal gas at T = Tbec) vs. T/Tbec- The value of 7 cr has been increased to 
more clearly reveal the anti-crossing behavior at low temperatures. As discussed in 
Section IV, the low temperature results only indicate the qualitative behavior. 
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